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MULTICARMER MODULATION USING WEIGHTED PROTOTYPE 



1 - Field of the Invention 

1-1 General Points 

The field of the invention is that of the transmission or broadcasting of digital 
data, or of analog and sampled data, designed to be received especially by mobiles. 
More specifically, the invention relates to the implementation of OFDM/OQAM type 
multicarrier signals. In other words, the invention applies to density 2 or even higher 
density signals. 

It is known that multicarrier modulation has many useful features, especially 
when it is associated with error-correcting encoding and interleaving. The COFDM 
(Coded Orthogonal Frequency Division Multiplexing) technique has also been chosen 
for the European digital audio broadcasting (DAB) standard and for the terrestrial 
digital video broadcasting (DVB-T). 

The COFDM technique offers a particularly simple system of equalization, 
namely the use of a guard interval, also called a cyclical prefix. This cyclical prefix 
provides for robust behaviour in the face of the echoes but at the cost of a relatively 
major loss of spectral efficiency. 

This problem is discussed inter alia in the French patent application No. FR- 
95 05455 (in which the COFDM modulation is called an OFDM/QAM modulation). 
To overcome this problem, this patent document presents a new technique for the 
implementation of OFDM/OQAM type multicarrier modulations. 

It will be noted that the different types of modulation discussed hereinafter are 
designated in a slightly different way in this prior art document and in the present 
patent appUcation. The following table gives the correspondence : 
FR95 05455 : Present Document: 



FUNCTIONS 



OFDM/QAM 
OFDM/OQAM 



OFDM/QAM/OFDM 
OFDM/OQAM/NYQUIST 
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OFDM/OMSK 



OFDM/OQAM/MSK 
OFDM/OQAM/IOTA 



OFDM/IOTA 



The term terme « OQAM » refers to the « Offset Quadratic Amplitude 
Modulation » definition. This expresses the fact that, for such modulations, there is a 
temporal offset between the transmission of the in-phase part and that of the in- 
quadrature part of a QAM symbol. 
5 1-2 Applications 

The invention can be applied in very many fields, especially when high 
spectral efficiency is sought and when the channel is highly non-stationary. 

A first category of applications relates to terrestrial digital radio-broadcasting, 
whether of images, sound and/or data. In particular, the invention can be applied to 
10 synchronous broadcasting which intrinsically generates long-term multiple paths. It 
can also be advantageously applied to broadcasting towards mobile bodies. 

Another category of applications relates to digital radiocommunications. The 
invention can be applied especially in systems of digital communications towards 
mobiles using high bit rates. 
15 2 - Reminders 

2-1 Transmission channel 

In a radiomobile environment, the transmitted wave undergoes multiple 
reflections, and the receiver therefore receives a sum of versions delayed by the 
transmitted signal. Each of these versions is attenuated and phase-shifted randomly. 
20 Since the receiver is assumed to be in motion, the Doppler effect acts also on each 
path. 

The conjunction of these efforts results in a non-stationary channel with deep 
fading at certain frequencies (frequency selective channel). For the applications 
referred to here, the transmission band is greater than the coherence band of the 
25 channel (the band for which the frequency response to the channel may be considered 
to be constant on a given duration), and fading thus appears in the band, i.e. at a given 



point in time, certain frequencies are highly attenuated. 

2-2 Description of a multicarrier modulation 

A multicarrier modulation is above all a digital modulation, namely a method 

for the generation of an electromagnetic signal out of the digital information to be 
5 transmitted. The originality and value of such a modulation is that it subdivides the 

limited band allocated to the signal into subbands. In these subbands, which have a 

chosen width smaller than the coherence band of the channel, the channel may be 

considered to be constant for a duration of transmission of a symbol, chosen to be 

smaller than the coherence time of the channel. 
Q 10 The digital information to be transmitted during this period is then distributed 

^ over each of these subbands. This has two uses in particular: 

J - reducing the modulation speed (namely increasing the symbol duration) 

=fi without modifying the transmitted bit rate, 

zl - simply modelling the action of the channel on each of the subbands: 

s 15 complex multiplier. 

^ It will be noted that, in reception, a system of low complexity for the 

^ correction of the data elements received (complex division by the estimated channel) 

p enables a recovery of the information transmitted on each of the carriers satisfactorily 

except for the carriers that have undergone a deep fading. In this case, if no steps for 
20 protecting the information have been taken, the data elements conveyed by these 
carriers will be lost. A multicarrier system therefore ensures that the generation of the 
electrical system must be preceded by digital data processing (error corrective 
encoding and interleaving). 

The patent No. FR 95/05455 gives a detailed description of the two types of 
25 existing multicarrier modulation. Their characteristics may be briefly recalled here. 
2-2-2 Notations 

Spacing between two adjacent carriers of the multiplex of carriers: Vq. 
Temporal spacing between two muhicarrier symbols transmitted (symbol 
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time): Tq. 
2-2-3 The prototype function 

The shaping filter for each of the carriers of the muUiplex is the same. It 
corresponds to the prototype function characterizing the multicarrier modulation. 

Let g(t) be this prototype function, the signal transmitted at each instant utq, 
on the m*^ central frequency subband v^, is a^^^^^e^"''" e'^^^"'^g(t - hIq) . 

In baseband, the expression of the signal transmitted, centered on the 
frequency MVq is therefore: 
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n m=0 

The functions e*^'"'" e^^'^'^^^^git- hTq) are called time-frequency translated 
functions of g(t). To retrieve the information transmitted by each of the subcarriers, it 
is necessary to choose g(t) so that its time-frequency translated functions are 
15 separable. To be sure of this, it is laid down that these translated functions should be 
orthogonal in the sense of a scalar product defined on all the finite energy fiinctions 
(which is a Hilbert space in the mathematical sense). This space accepts two possible 
scalar products, namely: 

- the complex PS = \x{t)y{t)dt 

R 

f \ 

20 - the real PS \x\y) = \x(t)y{t)dt 

V/? J 

Thus two types of multicarrier modulation are defined: 

- complex type, or again OFDM/QAM: the function g(t) chosen guarantees an 
orthogonality of its translated functions in the complex sense (example: 
OFDM, also called OFDM/QAM/OFDM). In this case, (pm^ri = 0 md the data 

25 elements a^^ ^ are complex, 

- real type or again OFDM/OQAM: the function g(t) chosen guarantees an 
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orthogonality of its translated functions in the real sense (examples: 

OFDM/OQAM/NYQUIST, OFDM/OQAM/MSK, OFDM/OQAM/IOTA). 

In this case, (pm^n ^ (n/2)'^(m^n) and the data elements ^ are real 
2-2-4 Density of the "time-frequency" network 
5 Since these modulations are designed for transmission at high bit rates, they 

will be associated with fairly high spectral efficiency values, in the range of 
4bits/s/Hz (digital television). The mapping of the bits coming from the error 
correction encoder will thus be of the QAM type. 

For an OFDM/QAM modulation, the real and imaginary parts of a complex 
10 ftinction derived from the QAM constellation are transmitted simultaneously at every 
symbol time T^. 

In the case of an OFDM/OQAM type modulation, they are transmitted with a 
temporal offset (QAM offset) of half a symbol time (T3/2). For one and the same 
transmission band and one and the same number of subcarriers, it is therefore 
15 necessary, for one and the same bit rate to be transmitted, that the rate of transmission 
of OFDM/OQAM type multicarrier symbols should be twice that of the OFDM/QAM 
type multicarrier symbols. 

These two modes of transmission of information are characterized by the 
density of the time-frequency network d = l/(Vo Xo). The OFDM/OQAM type 
20 multicarrier modulations the correspond to a density d = 2, and those of the 
OFDM/QAM type correspond to a density of d = 1 . 

It may be noted that a multicarrier modulation is characterized by: 

- the density of the "time- frequency" network on which it is defined, 

- the prototype fiinction. 

25 The implementation of an OFDM/OQAM modulation with a density 2, and of 

the corresponding demodulation, requires substantial computation power and high 
storage capacity. This therefore underlines the fact that the corresponding 
instruments are relatively complex and costly. 



3 - Goals of the invention 

The invention is designed especially to overcome or reduce these drawbacks. 
More specifically, a goal of the invention is to provide techniques for the modulation 
and demodulation of the multicarrier signals that are simple and cost little to 
5 implement as compared with the direct approaches. 

In other words, it is a goal of the invention to give modulation and 
demodulation techniques of this kind that limit the number of operations to be 
performed and the necessary storage capacity. 

4 - Description of the invention 

10 These goals as well as others that shall appear hereinafter are achieved 

according to the invention by means of a method for the modulation of a multicarrier 
signal with a density 1/(vo.to)=2, formed by successive symbols, each comprising M 
samples to be transmitted and constituted by a set of 2M orthogonal carrier 
frequencies in the real sense, 

15 the interval between two neighboring carrier frequencies being equal to vq and the 
interval between the times of transmission of two consecutive symbols, or the symbol 
time, being equal to to, 

each of said carrier frequencies being modulated according to one and the same 
modulation prototype fiinction g(t) with a truncation length of 2Lto, 
20 the method comprising, at each symbol time, the following steps: 

- the obtaining of a set of 2M complex coefficients representing data to be 
transmitted; 

- the computing of 2LM linear combinations from said 2M complex 
coefficients obtained, the weighting coefficients used in these combinations 

25 representing said prototype fiinction g(t), so as to obtain 2LM coefficients; 

- the summing of said 2LM coefficients weighted in the predetermined storage 
locations of a memory comprising 2LM storage locations representing 2L 
groups of M distinct partial sums, 
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so as to gradually form, in said storage locations, over a duration of 2Lto, M 

samples to be transmitted; 

- the sending of said samples to be transmitted. 
Thus, according to the invention, the data elements to be processed are stored after 
5 weighting and not before. It is thus possible to reduce the memory capacity needed as 
well as the number of operations performed. The samples to be transmitted are built 
gradually in each field of storage. 

According to one advantageous embodiment of the invention, a sample to be 
transmitted at the instant Jiq + kxo/M, referenced s^+jm is written as follows: 

10 

21-1 

^k^-jM = X ^k,qQj-q -^^k,qQ+MJ-q ] 
q^O 



where: 

Cqj to C2M.1J are the 2M complex coefficients generated between the instants 
15 jxo and G+IK; 

q and pk,q said weighting coefficients. 
In the case of an OFDM/OQAM modulation, we will generally have: 
ocicq = 0 for q as an odd parity number; 
Pk,q = 0 for q as an even parity number. 
20 The number of operations performed is therefore further reduced by half. 

In a preferred embodiment of the invention, the method comprises, for the 
generation of a symbol with an index j formed by M samples, the following steps: 

- the obtaining of 2M real inputs a^^ representing a source signal; 

- the pre-modulation of each of said real inputs producing 2M complex 
25 coefficients; 

- the reverse Fourier transform of said 2M complex coefficients producing 2M 
complex transform coefficients Qj to C2M-1 

- for each of the M pairs (C^j, C(k+M)j) of said transform coefficients, the 




8 




computation of 2L weighted coefficients, the weighing coefficients 
representing said prototype function; 

- the addition of the result of each of said weighted 2LM values to the contents 
of the 2LM distinct memory zones so as to gradually build the samples to be 
transmitted constituting the symbols j, (j+l), (j+2),... (j+2L-l); 

- the sending of M samples corresponding to the M oldest contents of said 
memory zones and then the resetting of the contents of said M memory zones. 
In general, said steps will be implemented at the rate Tq/M of the samples. 
The checking of the storage means is very simple. Thus, each sending step 

may be followed by a step for the updating of said memory locations comprising: 

- a physical shifting of the contents of each of said memory locations if the 
latter are elements of a shift register; or 

- an updating of the write and read addresses of said memory locations, if the 
latter are elements of a RAM. 

According to an advantageous characteristic of the invention, said coefficients 
representing data elements to be transmitted are obtained by the implementation of a 
mathematical transform comprising the following steps: 

- the apphcation of a real reverse Fourier transform; 

- the circular permutation of the result of this reverse transform by M/2 
coefficients leftwards; 

- the multiplication of each of said coefficients by i". 
It is thus possible to obtain complex transform coefficients from an FFT with 
real inputs. Again, this makes it possible to limit the number of operations 
performed. 

It is fiirthermore possible to simplify the computations by sUghtly modifying 
the equation of the signal centered on the frequency Mvq so that it is written as 
follows: 
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The invention also relates to the modulation devices implementing a 
modulation method of this kind. 

According to a particular embodiment, this device comprises especially: 

- means of mathematical transformation delivering said coefficients 
5 representing data elements to be transmitted at the rate xJTM. and in the 

following order (Coj, Cm+ij)v..(Cm-ij, Qm-ij); 

- 2LM-M simultaneous read/write RAM type memory locations; 

- N complex multipliers working at the rate Nto/2LM, N being equal to 1,2, 
4,..,or 2L. 

10 Thus, the memory space is further reduced. 

The invention also relates to a method for the demodulation of a received 
signal corresponding to a multicarrier emitted signal with a density 1/(vo.to) = 2. 
According to this method, an estimation of 2M real data elements transmitted at a 
given symbol time is reconstituted by means of the following steps: 
15 - the sampling of said signal received at the sample frequency Tq/M, delivering 

M complex samples received; 

- the storage of each of said M complex samples received in a predetermined 
location of an input memory comprising 2ML complex locations, in which 
there have been previously memorized (2L-1)M samples received during the 

20 21-1 previous symbol times; 

- the multiplication of the 2ML values contained in said input memory by 
coefficients representing said prototype function; 

- temporal aliasing, by the summing up of 2M series of L results of 
multiplication, so as to obtain 2M complex values; 

25 - the processing of said 2M complex values to form said estimations of the 2M 

real data elements transmitted. 

Advantageously, the 2M complex values derived from the temporal aliasing 
step between the instants (j+2L-l)To and (j+2L)To are written as follows: 



2L-1 

<t =0 
2Z,-1 

q^=0 

where: 

5 I'k'+o+q')!^ represents the sample received at the instant k'Xo+(j'^q')T^{/M; 

a\ q and P'^^, are said weighting coefficients. 
Most usually, the computations will be simplified because: 
^ k,q' ^ 0 for q' as an odd parity value; 
P'm' ^ 0 for q* as an even parity value. 
10 According to a preferred embodiment, said processing step comprises the 

following steps: 

- the application of a mathematical transformation that is the reverse of the 
one performed during the modulation on said 2M complex values delivering 
2M transformed values; 

15 - the correction of phase and/or amplitude distortions due to the transmission 

channel; 

- the extraction of the real part of said transformed complex values. 
In general, said steps are implemented at the rate x/M of the samples. 

The invention also relates to the demodulation devices implementing this 
20 method. These devices comprise: 

- means for the sampHng of said received signal; 

- means for the temporary storage of the complex sample functions 
comprising 2ML complex locations; 

- means for the multiplication of said memorized samples by weighting 
25 coefficients representing said prototype function; 

- temporal aliasing means summing up L weighting results so as to obtain 2M 



complex values; 

- means for the processing of said complex values delivering an estimation of 
2M real data elements transmitted at each symbol time. 

It is possible to further reduce the memory space needed in this device by 
5 means of: 

- storage means comprising 2ML-M simultaneous write/read RAM type 
complex memory locations; 

- N complex multipliers working at the Nto/2LM rate, where N is equal to 1, 
2,4... or 2L; 

10 - means of mathematical transformation working at the xJlM rate, whose 

inputs R^jto R2M.1J are read in the order (Roj, Rmj), (R,j, Rm^ij).-. (Rm-ij. R2M- 
io)- 

Finally, the invention more generally relates to a filtering method delivering 
series of M complex output values at regular intervals fi*om 2L series of 2M complex 
15 input values, 

said M complex values corresponding to a weighted sum of 2L of said complex input 
values to be processed, 

comprising the following steps for each series of complex input values: 

- the computation of 2LM linear combinations from said 2M complex 
20 coefficients obtained, the weighting coefficients being derived fi-om 2L real or 

complex filters with a size M, 
so as to obtain 2LM coefficients; 

- the summing of each of the weighted values in a predetermined memory 
location out of a set of 2ML memory locations each containing a partial sum 

25 so as to gradually form said output values in said memory locations on a 

period corresponding to the reception of 2L series of complex input values. 
The term "filtering" must of course be taken here in its general sense of 
processing or computation performed on data elements. This processing which 



i i 



comprises the computation of a weighted sum is done gradually. 
5 - Description of a preferred embodiment 
5-1 List of figtires 

Other features and advantages of the invention shall appear more clearly from 
5 the following aescription of a preferred embodiment of the invention given as a 
simple non-restrifitive illustration, and from the appended drawings, of which: 

Figure 1 gives a general and simplified illustration of the method of 
modulation ©f the invention (step j) used to generate M samples; 

Figure 2 illustrates the gradual construction of the M samples where 
1 0 Z=4; for the Io\a waveform; 

Figure 3 gives a more detailed illustration of the working of the 
method of modulation of the invention for the instants j-1 to j+2; 

Figure 4 Jmecifies the initiation of the procedure of modulation of 
Figure 3 where L 44; 

15 - Figure 5 is ay schematic diagram of a complex IFFT circuit known per 

se; 

Figures 6A to\6C illustrate the optimized architectures implementing 
FIFO systems and respectively using a single multiplier (Figure 6A), L 
multipliers (Figure 6B) or 2L multipliers (Figure 6C); 
20 - Figure 7 shows ai\ optimized embodiment of the reverse FFT using a 

real input FFT; 

Figure 8 illustrates \he working of the demodulation method of the 
invention when Z=4; 

Figure 9 shows the general case of demodulation deduced directly 
25 from Figure 8; 

Figure 10 illustrates a coiVesponding demodulator architecture; 

Figures 11 to 12 show two modes of implementation of the reception 
filtering in the case of a FIFO\ structure respectively using L and 2L 



i 



13 




multif^liers. 
5-2 Notations 

. Intercarrier spacing: Vq. 

. Intersymbol duration: Tq. 

. Density of the network: 1/(VoTo) = 2. 

. Band allocated to the signal: W= 2Mvq, 

. SampHng frequency: = l/T^ ^M/Tq. 

. Length of truncation of the protocol function: 2Ltq, 

In theory, the prototype function is the temporal support and/or infinite 
sequential support. However, to implement the corresponding digital filter, 
this function must be truncated. 

This is the case in OFDM/OQAM/NYQUIST (infinite temporal 
support) and in OFDM/OQAM/IOTA (infinite temporal and frequency 
supports). Typically, for the function Iota, L = 4 at the minimum. For a 
sampling at T^ defined here above, the digital function will have the length of 
2LM real coefficients. 
. Indices of the samples: 

To be consistent with the formula given hereinafter, the following is 



the prototype function the emitted signal s(t) and the emitted signal r(t) 

may be substituted foTx(t). 
5-3 Modulation algorithm 
5-3-1 Principle 

The baseband signal centered on the frequency Mv^ = fJ2, is written as 
follows: 



noted: 




(2) 
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^(t) = ZuLOn,J e "git -nr.) 

n m=0 (j) 



A sample of the signal is therefore written as: 

^ « m=0 ^ 



With the notations introduced here above (formula (2)), we have, after 
computation: 



2M-\ 2i^^ 



n m =0 



Given the rapid decrease of the prototype function, only the samples indexed 0 
to 2ML-1 are considered to be non-zero. We should therefore have 
0 < p- nM < 2ML - 1 . Taking p = k-h jM where 0 <k<M-\, we obtain 
15 j-{2L-\) < n <j\ 

The equation (5) becomes: 



2. ^m^n' (-1) ^ gk^U-n)M (6) 



/7= J~{2L-\) m = 0 



20 Finally, assuming q = j-n,wQ obtain the formula from which we derive the 

modulation algorithm: 
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2L-1 2M-1 

n n 1 ^'4-^4^4 4 4^2 4 4 4 4 4 4 4^ 



m (k + q M) 

liTt — ^ - 

2M 



with 



(7) 



14 4 4 4 4 4 4 ^J^^^f^^WZ^ 444444444444^ 
144444444444 4Rfv^r4e4FfT4 44444444444444 43 
Weighting by the prototype function g(t) 

0<A:<M-1 



Noting c^nj-q pre-modulated input values giving: 



-m+ij-q). ^>jn(j-q+L) 
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the equation (7) becomes: 



2L-\ 
q = 0 
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X ^fn,j-q ^ 



m {k+q M)^ 



2M 



m = 0 



gk+qM (8) 



This formula leads to a modulation algorithm in three main steps: 

- Pre-modulation of the data elements by means of a simple complex 
15 multiplication. 

- Reverse Fourier transform (by IFFT algorithm). 

- Filtering by the prototype function. 

Here, as in the rest of the document, the term filtering is understood to mean 
an operation of weighting of the results of 2L reverse FFT operations by 



certain values of the prototype function, followed by an operation of 
summing of these weighted coefficients. In other words, this is a linear 
combination. 

Here below, we present two possible modes of implementation. The second 
5 mode is optimal and is more precisely the object of the invention. 

Although in practice, the work is done at the sample rate t^/M, we shall retain 
the block structure of M samples to describe these modes of implementation with 
greater clarity. 
5-3-2 Direct architecture 
10 In view of the formula (8), it is necessary to perform IL complex IFFT 

operations with a size IM to generate M samples (corresponding to the duration of a 
multicarrier symbol Tq). 

However, the result of an IFFT comes into play on the computation of IL 
consecutive blocks of samples. To compute the M current samples, it is therefore 
15 necessary to compute only the reverse Fourier transform of the IM last data elements 
entered into the modulator, the results of the 21-1 other IFFT values having been 
computed in the previous steps and stored in the memory. 

The modulation algorithm therefore comprises the following main steps: 

- The pre-modulation of the 2Mreal inputs delivering 2M complex values. 

20 - The reverse Fourier transform with a size of 2M complex values (IFFT 

algorithm). 

- Storage of the result by re-updating a buffer with a size of 2L*2M complex 
values (containing the results of the IL reverse FFT operations indicated in the 
computation). 

25 - Filtering of the ILM elements of the storage buffer by the prototype function. 

- Sending of the M complex samples thus computed. 
The requisite memory size is therefore: 

- a RAM with a size of IL^IM complex values (input buffer), 



- a ROM with a size of real values (weighting coefficients). 

(The prototype function is chosen to be symmetrical, on the 2ML 
weighting coefficients, and only the ML+\ values are distinct.) 
This first procedure reveals a waste of RAM type memory. The second 
5 architecture proposed shows that it is possible to reduce the size of the necessary 
RAM by more than half This reduction is accompanied by a reduction in the number 
of operations and therefore an increase in the processing speed. 
5-3-3 Optimized architecture 

It is possible to optimize the architecture of the modulator according to the invention, 
10 both at the level of the filtering by the prototype function and that of the reverse FFT. 
Indeed, in analyzing the formula (8), it can be ascertained that for each of the 2L 
IFFTs involved in the computation of the current block of M samples, only M points 
on 2M are used. 

It is then possible to reduce the required RAM type memory by half by storing 
15 the data elements used for the computations of the different blocks of samples after 
filtering, and not before. Furthermore, the specific structure of the complex data 
elements at input of the reverse FFT (a^^ ^.|^ni(n+L) ^m+n^ enables the use of a reverse 
FFT algorithm with real inputs. 

In order to specify this method, we shall develop the formula (8): 

20 

^k+JM- I 2-f ^m,j-2q'^ Sk+lq' M 2^ ^ m, j-(2q'-i-\) ^ ^Jt+(2 ' +1 ) A/ I 

q' = 0 \_m=0 m=0 J 

Let: 

Q^-Z^.,.^ (9-1) 

m = 0 

25 and 

2M-\ 2i.T^-^^i^ 

C,,M^ = Z c„^„e ' 2^ (9.2) 

m = 0 



We have: 



^k+JM ~ ^ \f^kJ-2q'Sk+2q'M "^^k+M J-{2q'+\)S k+{2q'+\)M\ /i ^.x 
q'=0 \^^) 

with 0<k<M'-l 

5 The equation (10) expresses the construction of M complex values from 2ML 

complex values. It can be written more generally: 

2L-1 

^ = 0 ^ ^ 

with 0<k< M-l 

10 In the embodiment described, 

10 if q is odd - parity [sk +qM ^ odd - parity 

• o - and Bv n = ^ 

gk+qM if q IS even-panty ^^'4 q -f ^ gven-parity 

To generate M samples according to this modulation algorithm, it is therefore 
possible to proceed as illustrated in Figure 1 (step j): 
1 5 - Pre-modulation 1 1 of the 2M real inputs. 

Reverse Fourier transform 12 of the 2M complex data elements thus 
obtained so as to generate and C^^^^j . 

A weighting 13 (corresponding to the application of the prototype 
fimction) of the result of the reverse Fourier transform by the prototype 
20 function: L parallel weighting operations. 

The L weighting vectors, with a size 2M, have the following coefficients: 

[gov, giM-xl [gzMv, g4M-i]v..,[g2i:M-2M,", gzLM-i]- 
The addition 14 of these weighting results to the output buffer of with 
a size of 2ML.complex values 
25 - The shifting 15 of the output buffer with the sending of M samples, 



corresponding to the M oldest values contained in the buffer. 

A sample of the signal to be transmitted represents a sum of 2L weighted IFFT 
results. Each block of M consecutive memory slots of the output buffer contains M 
partial sums of 2L-m^^^^^ terms each, where m^iock varies from 1 to 2L (m^iock =2Z 
5 corresponds to the "all at zero" block, due to the buffer shift operation (step 14) at the 
instant 0-1)To). The 2ML elements coming from the L parallel weighting operations 
are herein added to the 2ML elements of the buffer. 

After this operation, the M partial sums of the buffer corresponding to m^iock = 
1 are completed, namely the M current samples are computed and may therefore be 
10 transmitted. 

This operation is described in Figure 2, where Z=4. Each line illustrates the 
situation of the construction buffer of the data elements to be transmitted, at a given 
instant. It is necessary to have 2L consecutive symbol times to gradually construct a 
sample to be transmitted. 
15 The waveform 2L shown corresponds to the Iota fianction. It is represented by 

the 2L vectors of coefficients 22 [gj to [g^Ml^ where the index k varies from 0 to M- 
1. 

At each instant, the 2M coefficients at input are multiplied (23) by the 
coefficients 22 and then added up (24) each to a partial sum. 
20 The M partial sums complemented at the step 15 are transmitted, the contents 

of the buffer are shifted by M memory slots (so as to ensure the right order of 
computation of the next M samples) and M zeros are inserted in the M vacant memory 
slots. 

The diagonal 25 thus illustrates the computation of Sj,+ji^. 
25 Figure 3 gives a more detailed view of the working of this algorithm for the 

instants j-1 to j+2. If we consider the instant j, the coefficient a^ to be transmitted 
suppUes the pre-modulation module 31, which gives the reverse FFT 32 the 
coefficients c^g. The reverse FFT delivers the and C^^^^ values (the index k varies 
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from 0 to M-1) subjected to the weighting 33 (the weighting operations in parallel) to 
deliver the results 34 which are summed up in an output buffer 35. 

Figure 3 gives an indication of the exact contents of these output buffers. 

Figure 4 illustrates the triggering of this modulation procedure where L = 4. 
5 The architecture of the modulator corresponding to the above algorithm 

presented here above must therefore comprise: 

a ROM with a size of real values containing the coefficients of the 

filter, 

a RAM with a size oilML complex values corresponding to the output buffer, 
10 a complex FFT circuit (achieving a reverse FFT) with a size 2M. 

To increase the processing speed, the weighting operations will be made 
parallel by using L RAMs with a size 2M associated with L multipliers or even 2L 
RAMs with a size of M complex values associated with 2L multipliers instead of one 
RAM with 2ML complex values. 
1 5 The complexity of the modulation circuit is therefore: 

for the filtering: 

In order to carry out the filtering, we multiply the results of the complex 
reverse FFT by the 2ML coefficients of the prototype function, by carrying out L 
weighting operations in parallel, including the result of the output buffer. Given that 
20 the coefficients of the prototype function are real, we have (2x2ML) real 
multiplications and 2ML complex additions or AML real additions. The size of the 
output buffer is then 2ML complex memories or AML real memories. 

for the IFFT transform: 

The results indicated in the following table relate to a conventional complex 
25 BFFT circuit whose schematic diagram is given in Figure 5. It comprises an input 
buffer 51, receiving 2M inputs, a computation unit 52 supplied by coefficients stored 
in the ROM type memory 53 and computed values stored in a RAM type memory 54 
that deliver 2M outputs. A control module drives these different elements. 



This table gives an estimation of the operations and equipment needed for the 
modulation part: 



Modulation 


Addition 

Operations 

(real) 


Multiplication 

Operations 

(real) 


RAM 
(real) 


ROM 
(real) 


FIFO 
(real) 


Pre-modulation 


X 


X 


X 


X 


X 


Reverse FFT 


6M(l + log2 M) 


4M(l + log2 M) 


m 


2M 


X 


Filtering 


4ML 


AML 


AML 


ML+l 


X 



It will be noted that no multiplication or addition is needed for the pre- 



5 modulation because the simple complex multiplication to be found at this stage is 
expressed, at the level of the architecture, by permutations of real and imaginary parts 
as well as changes of sign. 

According to one variant of the invention, an additional gain in memory space 
may be obtained. It is indeed possible to use only (2L-1) RAMs with M complex 
10 values (total storage: 2ML-M complex values instead of 2A/L). To do this, at the step 
j, a reading is done of the M samples to be transmitted into the corresponding RAM 
gradually, and the M complex values Cf^^j^jgi^^2L-\)M written progressively at the 
same addresses. To carry out this filtering operation, it is possible, as required, to use 
RAMs or FIFOs. 

15 Figures 6A to 6C illustrate this method in the case where a FIFO structure is 

chosen. 

In the case of Figure 6 A, a single multiplier 61 is implemented. It multipUes 
data elements delivered by the reverse FFT by the weighting coefficients and supplies 
an adder 62 that also receives the output of the FIFO memory 63 containing 2ML-M 
20 complex values. This FIFO 63 is supplied by the result of the addition 62. A control 
module 64 enables the output of the FIFO to be directed outwards to deliver the M 
complex values ready to be transmitted 65. 



i / 



It is possible to use N parallel-connected multipliers, where N = 1,2,4,...2L. 

Thus, in the case of Figure 6B, L(=4) multipliers 61 j to 6I4 are implemented 
in parallel. They are supplied alternatively with one or the other of the weighting 
coefficients associated with them. 

Each of them supphes an adder 62, to 624 which also receives data elements 
from the 2L FIFO memories 63, to 63^ each comprising M complex values. The 
FIFO memory 63, delivers the M complex outputs. Selection means 66 enable the 
selection of a FIFO memory to be taken into account at each point in time. 

Figure 6C shows the implementation of 2L (=8) multipliers. In this case, the 
structure no longer requires the presence of control means. The 2L FIFO memories 
63, to 63^ are each supplied by a multiplier 61, to 6l8, associated with its own 
weighting coefficient and associated with an adder 62, to 625. 

It must be noted that the reduction of the memory space needed for the output 
buffer is valuable only if the following two conditions are met: 

- the algorithm used to obtain the reverse FFT arranges its output in the 
optimum order Coj,CMj,C,j,CM+,jv",CM.ij,C2M.i and works at the rate To/(2M). 
Indeed, in the case of outputs arranged in the reverse bit order or even in the 
natural order, Coj,C,j,...C2M-2j5C2M-i the proposed gain in memory due to this 
simultaneous reading and writing of the output buffer then requires the 
reordering of the outputs of the reverse FFT in the optimum order which 
requires a storage at output of the FFT. In either case, the gain in memory 
will be negligible or even zero. 

- The multipliers used, to prevent any storage of the outputs of the reverse 
FFT, work at a high speed: the rate of operation of N parallel-connected 
multipliers must be equal to Nto/(2LM), for N = 1,2,4,...2L. 

According to another variant of the invention, it is possible to optimize the 
revere FFT. Given the « particular » complex character of the inputs (a^n(-l)"^"^^^i"^") 
of the reverse FFT at transmission, it is possible to use an FFT algorithm with real 



inputs. 

It is known that phase-shifting the inputs of an FFT with a size 2M, by i"" amounts 
to applying a circular permutation to its outputs y^ by Mil leftwards. By applying 
this result, it can be clearly seen that the pre-modulation step: 

followed by the complex reverse FFT can be done as illustrated in Figure 7 in the 
following steps: 

- real reverse FFT 71 of the data elements (.1)"^^"^^) ^ ^ 

- circular permutation 72 by Mil of the outputs, 
10 - multiplication 73 by i". 

An algorithm of this kind enables a reduction by half of the memory space 
needed for the FFT as well as the number of operations. Figure 5 shows these three 
operations. 

It can be noted that the operation of multipUcation by (-i)'"^"^^) has been 
15 omitted. Indeed, it can be avoided. 

At transmission, to generate the baseband signal s(t) (equation 1) digitally, the 
values of a^^ must be multiplied by in reception, the estimation of the data 

requires this multiplication again as shall be seen hereinafter. 

Given the fact that the withdrawal of this multiplier term has no effect on the 
20 orthogonality of the time-frequency translated values, it is possible to remove the 
need for this multiplication. This amounts then to generating the baseband signal 
centered on the frequency Mvq = ijl according to: 



25 



2A/-1 
n m=0 



6 - Demodulation 

The method of demodulation must enable a recovery of the useful information 
transmitted through the samples of the signal received in reception. It is assumed 



here that the "Doppler-delay" channel (the most general case) of the transfer fUnction 
T(f,t) is perfectly estimated and that it is locally likened to a complex multiplier 
channel T„,j = p^je^""' • 

Given the orthogonality of the "time-frequency" translated functions of the 
prototype function, the information sent at the instant Jiq, on the carrier m is thus 
estimated: 



^.j = e"''-^- \r(t)gljm (12) 

\_PmJ 



Tr) 1 

10 In practice, we work on the versions sampled at — = — of the received 

M W 

signal, the demodulation function then becomes: 



W 



<(P-yAo-^J (13) 



1 



15 Resuming the notations given by the formula (2) and taking account of the 

limited number of coefficients representing the prototype function {2ML), we obtain a 
demodulation formula as follows: 



m,j 



^ jM+2gM Sk+2qM 

Ij't 4 42 4 4 43 



-2ifun- 



2M 



with 



1 (_l)'«0>^) 

o W 
f'^- 4 4 41 2 4 4 4 4 « 

Phase and amplitude correction | W^ig)jti^q ^Ipy ^^h^^^tppypje ^n^ti^n /j 4^ 

Complex FFT 

p = jM + IqM + k 

0<k<2M~\ 

Q<q<L-l 



The formula (14) suggests five steps for the fast demodulation algorithm.: 

- weighting of samples received by the prototype function, 

- temporal aliasing, 

- 2M sized complex FTT, 

- phase and amplitude correction, 

- extraction of the real part. 

The solution proposed according to the invention to minimize the memory 
space taken up by the filtering (weighting and temporal aliasing) in reception is as 
follows: 

- inserting the M samples received in an input buffer with a size of 2ML 
complex FFT functions, 

- multiplying the data elements of this buffer by the coefficients representing 
the prototype function, 

- summing the results of these multiplications (temporal aliasing), 

- applying a direct Fourier transform to the 2M complex values thus obtained, 

- correcting the result of this FFT in phase and amplitude, 

- extracting the real part. 

It must be noted that this technique is independent of the way in which the 
signal has been constructed at transmission. It can be applied to the reception of any 
type of OFDM/OQAM multicarrier signal. 

In order to illustrate the working of this algorithm, we shall break down the 
complex 2M sized complex signals into two M sized complex sub-signals each, as 
follows: 

^k+JM+2qM ~ 'jt'+yA/+2i7A/ ^ ^{)t£{0,...^-l}} ^{k'+M)+JM+2qM ^^{ke{M ,...,2M-\}} 




(15.1) 



and 



Sk+jM+2qM Sk'+jM+2qM ^ l{jt€{o,...^-l}} S(k'+M)+jM+2qM ^ ^{ke{Af ,.,.,2M -\}] 

\0<k<2M-l (15.2) 
with < 

\0<k'<M-l 
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The entry of FFT for its part will be referenced \j = Ya^k^jM^iqM^gk^iqM ^ k 
going from 0 to 2M-1. 

5 The estimation of the data elements sent will start after a delay of (2Z-1)To. It 

is necessary indeed that all the received samples comprising the data element amj 
should be stored in the input buffer before dyyij, is computed. 

The notations (1 5. 1) and (15.2) here above make it possible, when Z=4 , the 
inputs of the FFT needed for the estimation of a^^o, to write these equations in the 
10 form (with k' = O...M-l): 



^k+M,0~ ^k' + M^k'+M'^ *^ 't'+3M^it'+3M ^' + 5A/?Jfc'+5 AY ^k'+lM^k'+lM 

Similarly, the inputs of the FFT corresponding to the estimation of the values 
15 a^ 1 (m=0 to 2M-1) are thus built: 

^k,\ ~ ^k'+M^k' ^k'+3M8k' + 2M 't' + 5A/&'+4M ~^ ^k' -hlMSk' +6 M 

^•k'+M,\~ ^k'+lM^k' + M ^k'+AM^k'+iM Ifc'+eA/^Jt' + SA/ ^k' +B M^k'+IM 

Figure 8 illustrates the working of the architecture proposed in the case where 
Z=4. The general case is illustrated by Figure 9 and is deduced directly from Figure 
20 8. 

In this figure 8, the M samples 81 received at an given point in time are stored 
in an input buffer 82. At each symbol time, the data elements contained in this buffer 
82 are multiplied (83) by the weighting coefficients 84 representing the waveform 85 
(IOTA in the example) then added up (86) to carry out the aliasing. 
25 The corresponding data elements R^^ supply the FFT 87, performed on 2M 



complex samples. Then a phase correction 88 is done and then an extraction 89 of the 
real part. 

Finally, a shift is made of the contents of the input buffer 82. Figure 8 
presents the contents of this buffer for eight successive instants, corresponding to the 
5 production of the outputs a^j a k^-^-j. 

The above formulae represent the case of the demodulator associated with an 
OFDM/OQAM modulator with a density 2. However, this architecture remains 
applicable to the case of the restitution of 2M complex values from M complex values 
resulting from the general case of modulation illustrated by the formula (11). The 
M 10 associated general formula would be: 

r~ 

£ r 2i-i 

S with A:'g{0,...,M-1} (16) 

ffl ^k'+MJ ~ ^^^k\J+q'-^ k\q' 

y The above algorithm requires the following means: 

p - a RAM with a size of 2ML complex values (input buffer), 

^ 15 - a ROM with a size of (ML+1) real values (coefficients of the digital 

a filter), 

- 1, L or 2L complex multipUers depending on the degree of parallel 
performance of the weighting operations, 
- a complex FFT circuit with a size 2M. 
20 Figure 10 illustrates the architecture proposed. 

The input buffer, capable of containing 2ML complex values, receives M 
samples at each tq. The weighting by the prototype ftinction is done by the 
multiplication operations 102 and then the aliasing is done by means of two adders 
103, and 1032, which supply a buffer 104 of 2M values supplying the complex FFT 
25 105. 

At the end of the FFT operation 105, a phase and amplitude correction 106 is 
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performed and then the real part is selected at 107 to give the 2M real values 
transmitted. 

Just as at transmission, it is possible to perform the weighting operations in 
parallel (complex multiplication) by using 2Z buffers with a size of M complex values 
5 associated with L, or 2L multipliers, rather than a single one with a size of 2ML. 
These aspects are illustrated respectively by Figures 11 and 12. The operation is 
deduced directly from the one described with reference to Figures 6B and 6C for 
transmission. 



10 2ML-M complex values instead of 2ML) to store the received samples. It is 

necessary, at the step j, to read the (k=O..M-l) sample received at the step (j"(2L- 

1)), and write the current sample at the same address. 

To perform this operation, it is possible, as needed, to use RAMs or FIFO 

memories. Figures 1 1 and 12 implement this method in the case of a FIFO structure. 
15 Once again, as at transmission, it must be noted that the reduction of the 

memory space needed for the input buffer is useful only if the following two 

conditions are met: 



It is possible to use only (2L-1) RAMs of M complex values (total storage: 



20 



- the algorithm used to achieve the FFT works at the rate xJ{2M) with the 
inputs arriving in the « optimum » order R^j,RMj,R,j,RM+ijv,RM-io.R2M-ij- If 
not, a storage at the input of the FFT will be necessary and the memory gain 
will then be negligible or even zero, 



- the multipliers used, also to prevent any storage of inputs of the FFT, work 
at high speed: the rate of operation of N muUipHers in parallel must be equal 
to Nto/(2LM), for N = 1,2,4,.. .2L. 
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